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ABSTRACT 

We study the acceleration of charged grains by niagnetohydrodynamics (MHD) turbulence in the 
interstellar medium (ISM) . We begin with revisiting gyroresonance acceleration by taking into account 
the fluctuations of grain guiding center along a uniform magnetic field (i.e. nonlinear theory-NLT). 
We calculate grain velocities due to gyroresonance by fast MHD modes using the NLT for different 
phases of the ISM, and compare with results obtained using quasi- linear theory (QLT). We find for the 
parameters applicable to the typical ISM phases that the fluctuations of grain guiding center reduce 
grain velocities by less than 15%, but they can be important for more special circumstances. We 
confirm that large grains can be accelerated to super- Alfvenic velocities through gyroresonance. For 
such super- Alfvenic grains, we investigate the effect of further acceleration via transit time damping 
(TTD) by fast modes. We find that due to the broadening of resonance condition in the NLT, the 
TTD acceleration is not only important for the cosines of grain pitch angle relative to the magnetic 
field fj. > Va/v, but also for fi < Va/v where v is the grain velocity and Va is the Alfven speed. 
We show that the TTD acceleration is dominant over the gyroresonance for large grains, and can 
increase substantially grain velocities induced by gyroresonance acceleration. We quantify another 
stochastic acceleration mechanism arising from low frequency Alfven waves. We discuss the range of 
applicability of the mechanisms and their implications. 

Subject headings: dust, extinction - ISM: kinematics and dynamics, acceleration - ISM 



I 1. INTRODUCTION 

Dust grains play crucial roles in many aspects of the in- 

] terstellar medium (ISM). For example, alignment of dust 
grains with respect to magnetic field provides insight into 
star formation through far-infrared and submm polarized 

[ emission (see Lazarian 2007 for a review). Very small 
spinning dust grains radiate microwave emission that 
contaminate to cosmic microwave background (CMB) 

I radiation (Draine & Lazarian 1998; Hoang, Drainc & 
Lazarian 2010; Hoang, Lazarian & Draine 2011). Opti- 
cal extinction and polarization properties of dust depend 
mainly on its size distribution. Grain-grain collisions, 
which depend on grain relative motions, govern grain co- 
agulation and destruction, that result in the grain size 
distribution (see Hirashita & Yan 2009). Grain-grain 
collisions arc considered the first stage of planetesimal 
formation in circumstellar disks (see e.g., Dullemond & 

' Dominik 2005). 

Traditionally it is believed that the motion of dust 
grains in the ISM arises from radiative force, ambipo- 
lar diffusion and hydrodrag (see Draine 2011). The re- 
sulting motion from these processes is sub- Alfvenic (i.e., 
V <^ Va = B/y^Airp where B is the magnetic field 
strength and p is the gas mass density) , except in special 
environment conditions (see Purcell 1969; Roberge et al. 
1993). 

Astrophysical environments are practically all magne- 
tized and turbulent, and turbulence is expected to be an 
important factor in accelerating dust grains. The evi- 
dence for turbulence from electron density fluctuations 
testifies for the existence of the so-called Big Power Law 
in the Sky (Armstrong et al. 1995; Chepurnov & Lazar- 
ian 2010), while the fluctuations of velocity (see Lazar- 



ian 2009 and references therein) provide convincing ev- 
idences of the dynamic nature of the observed inhomo- 
geneities. Whether an environment is thermally domi- 
nated or magnetized dominated depends on the plasma 
f3 parameter, which is defined as the ratio of gas pressure 
to magnetic pressure 

^ STrnHfeBTgas ^ . . ( \ f Tg^\ f W fiG V 

^ ■ V30cm-3/ ^100 K; B J 

where nn is the gas density and Tgas is the gas temper- 
ature. 

Recent decade has been marked by substantial progress 
in understanding of MHD turbulence. This included gen- 
eralizing incompressible Alfvenic turbulence^ by Goldre- 
ich & Sridhar (1995) to realistically compressible media 
and successful testing of the compressible theory (Lith- 
wick & Goldreich 2001; Cho & Lazarian 2002, 2003; 
Kowal & Lazarian 2010). In what follows in describ- 
ing compressible MHD turbulence we shall be guided by 
the mode decomposition of MHD turbulence into Alfven, 
slow and, fast modes presented in Cho & Lazarian (2002, 
2003). 

Studies of grain acceleration for magnetized turbu- 
lent environments were initiated by Lazarian & Yan 
(2002) who dealt with the acceleration by incompressible 
Alfvenic turbulence. Comprehensive studies of the accel- 
eration in realistically compressible environments were 
performed in Yan & Lazarian (2003, hereafter YL03) and 

^ While there are still ongoing debates about the detailed struc- 
ture and dynamics of incompressible MHD turbulence, we believe 
that Goldreich & Sridhar (1995) model provides an adequate start- 
ing point. In fact, recent studies in Bercsnyak & Lazarian (2010) 
and Bercsnyak (2011) provided additional supports for the model. 
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Yan, Lazarian & Draine (2004, hereafter YLD04). Those 
studies identified gyroresonant interactions of grains with 
fast MHD modes as a new powerful mechanism of grain 
acceleration. Recently, Yan (2009) considered betatron 
acceleration and came to the conclusion that for most en- 
vironments the betatron acceleration is subdominant to 
the gyroresonance acceleration for sub-Alfvenic grains. 
The application of grain velocities predicted by gyrores- 
onance for modeling dust extinction curve provided good 
correspondences between observations and theoretical 
predictions (Hirashita & Yan 2009) indicating that the 
turbulence is indeed the main driving force behind grain 
acceleration. 

The gyroresonance acceleration due to compress- 
ible MHD turbulence (YL03;YLD04) was studied using 
quasi-linear theory (QLT, Jokipii 1966; Schlickeiser & 
Miller 1998). The underlying assumption of the QLT is 
that the guiding center is assumed to move in a regu- 
lar trajectory along a uniform magnetic field Bq. The 
condition for a grain with velocity v to resonantly in- 
teract with fast MHD modes at the scale k]\ is given by 
Lo — — uj = nil, for n = 0, ±1, ±2, ... where ^ is the 
cosine of the grain pitch angle between v and Bq, ui is 
the wave frequency, and 51 is the Larmor frequency of 
the charged grain around Bq. The gyroresonance accel- 
eration with n 7^ is dominant by eddies with size equal 
to gyro radius I rg. 

Transit-time damping (TTD) or transit-time accelera- 
tion, arises from resonant interactions of particles with 
the compressive component of magnetic fluctuations (i.e., 
the component parallel to the mean magnetic field Bq). 
When the grain moves together with the wave along Bq, 
it is subject to magnetic mirror forces —{mv\/2B)V\\B, 
where m is the grain mass, is the grain velocity com- 
ponent perpendicular to Bq, and B is the total magnetic 
field. In the plasma reference, the back and forth col- 
lisions of the grain with the moving magnetic mirrors 
increase grain energy because the head-on collisions are 
more frequent than trailing collisions due to the larger 
relative velocity between grain and wave (sec Fisk 1976; 
Schlickeiser & Miller 1998). The TTD acceleration with 
resonance condition fc||z;// = u, was disregarded in previ- 
ous studies on grain acceleration because grains are ex- 
pected to move slowly along the uniform magnetic field, 
for which they can not catch up with the propagation 
of magnetic mirrors along this direction. Although the 
gyroresonance is found to be able to accelerate large 
grains to super- Alfvenic velocities (see YL03; YLD04), 
which is sufficient to trigger TTD, the resulting grain 
motion mostly perpendicular to the uniform field Bq (i.e. 
u/i = 0) in the QLT regime makes TTD incapable. 

Due to magnetic fluctuations in the ISM, the local 
magnetic fleld B can be decomposed into a uniform fleld 
plus a turbulent component, i.e., B = Bg + (5B. Thus, 
any perturbation 5B will induce the fluctuations of grain 
guiding center from a regular trajectory along the uni- 
form magnetic field (see e.g., Shalchi 2005). Non-linear 
theory (hereafter NLT) for gyroresonance that takes into 
account such fiuctuations of guiding center was formu- 
lated in Yan & Lazarian (2008, hereafter YL08) to de- 
scribe the propagation of energetic particles, and later 
it was applied in Yan et al. (2008, hereafter YLP08) to 
study acceleration of energetic particles in solar fiares. 

The important modification present in the NLT is the 



broadening of resonance function from a Delta function 
(5(aj— fc||U|| — nl7) to a Gaussian function Rn{Lo~k\\v\\—n^^) 
(see YL08). Such a broadening of the resonance condi- 
tion allows grains moving with v Va perpendicular 
to Bq to have TTD with compressive waves propagating 
along Bq. We are going to clarify the effects that TTD 
induces on grain acceleration in the present paper. 

In what follows, we revisit the basics of resonance ac- 
celeration for charged grains by taking into account addi- 
tional physical processes that were not considered within 
original treatments. In discussing the gyroresonance ac- 
celeration, we are going to take into account the fluctua- 
tions of grain guiding center. In particular, we are going 
to investigate the efiiciency of TTD on grain acceleration 
in MHD turbulence. 

The structure of the paper is as follows. In §2, we 
present briefly the problem of grain charging, important 
dynamical timescales, and identify the range of grain size 
in which grain charge fluctuations are important. We 
revisit gyroresonance acceleration, and introduce TTD 
acceleration in §3. Grain velocities induced by gyrores- 
onance acceleration and TTD are presented in §4. §5 
is devoted for stochastic acceleration by low frequency 
Alfvcn waves. Discussion and summary are presented in 
§6 and 7, respectively. 

2. GRAIN CHARGING AND DYNAMICS 

2.1. Grain Charging 

Charging processes for a dust grain in the ISM consist 
of its sticking collisions with charged particles in plasma 
(Draine & Sutin 1985) and photo-emission induced by 
hv > 13.6 eV photons (Weingartner & Draine 2001). 
In the former case, the grain acquires charge by cap- 
turing electrons and ions from the plasma, while in the 
latter case the grain looses charge by emitting photoelec- 
trons. After a sufficient time, these processes result in a 
statistical equilibrium of ionization, and the grain has 
a mean charge, denoted by (Z), which is equal to the 
charge averaged over time. Due to the discrete nature 
of charging events, the grain charge fluctuates around 
(Z). The probability of finding the grain with charge Ze 
is described by charge distribution function fz ■ Here we 
find the charge distribution fz using statistical ionization 
equilibrium as in Draine & Sutin (1985) and Weingartner 
& Draine (2001). Hoang & Lazarian (2011) found that 
the statistical ionization equilibrium is not applicable for 
tiny grains with size a < 10 A for which the charging is 
infrequent, but it is adequate for grains considered in the 
present paper. 

Figure 2 shows the variation of the grain mean charge 
\{Z)\ for graphite and silicate grains in the cold neutral 
medium (CNM), warm neutral medium (WNM), warm 
ionized medium (WIM). In the WIM, (Z) varies rapidly 
with the grain size, and change its sign at a ^ 10~^ and 
^ 10~^ cm. marked by filled circles. 

Let us define a characteristic relaxation time of the 
charge fluctuations, tz, which is equal to the time re- 
quired for the grain charge to relax from Z to the equi- 
librium state (Draine & Lazarian 1998b): 

^ {{Z-jZ))') ^ al 

J2z fzJtot{Z) J2z fzJtotiZ)' 

where Jtot {Z) is the total charging rate due to collisional 
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Fig. 1. — Mean grain charge \{Z)\ as functions of grain size a for 
graphite and sihcate grains in different ISM phases. For the WIM, 
(Z) I changes rapidly with a, and filled circles mark the change in 
grain charge between being positively charged (+) and negatively 
charged (-). 

charging and photocmission (see Draine & Sutin 1987; 
Weingartncr & Drainc 2001). Here we averaged over all 
possible charge states Z to find tz- 

In an ambient magnetic field B, the grain with mean 
charge {Z)e gyrates about B on a timescale equal to the 
Larmor period: 



TL 



2'Kmc 

\{zWb 



1.56 X 10^ 



10- 



;) {\{z)\b) 



yr, (3) 



where m = 4/37ra^pd with pd being the dust mass density 
is the grain mass. We adopt pd = 2.2 and 3.0 g cm~'^ for 
graphite and silicate grains, respectively. The Larmor 
frequency reads VL = {{Z)eB)/mc. 

We calculate the relaxation time of charge fluctuations 
Tz for both graphite and silicate grains in various phases 
of the ISM with physical parameters listed in Table 1. 
Figure 2 compares tz with the gas drag time Tdrag (see 
Eq. 5) and the Larmor period t^. It can be seen that 
Tz tl < Tdrag for grains larger than ~ 2 x 10"'' cm. 
For grains smaller than ^ 2 x 10~^ cm, tz > tl, so that 
the assumption for grains to have a constant charge is no 
longer valid. As a result, the fluctuations of grain charge 
should be accounted for in the treatment of resonance 
acceleration for such very small grains. This issue will be 
addressed in our future paper, in which we employ Monte 
Carlo method to simulate grain charge fluctuations (see 
e.g., Hoang & Lazarian 2011). In the present paper, for 
the sake of simplicity, we adopt {Z)e for grain charge 
within the entire range of the grain size distribution. 

2.2. Grain Translational Damping 

Interactions of dust grains with the ambient gas 
present the primary mechanism of dissipating transla- 
tional motions of grains. The damping rate of transla- 
tional motion arising from the interaction with neutral 
gas is essentially the inverse time for collisions with the 
mass of the gas equal that of a grain (Purcell 1969), 
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Fig. 2. — Timoscales for gas drag Tdrag, Larmor period r/,, and 
charge fluctuations tz as functions of the grain size a for subsonic 
silicate grains in the various phases of the ISM. Shaded area marks 
the range of grain size in which the charge fluctuations are impor- 
tant, i.e., Tz > Ti. The peaks in correspond to the change in 
sign of grain charge. 

where m„, n„, and T„ are the mass, volume density, and 
temperature of neutrals, and a is the grain radius. 

When the ionization degree is sufficiently high, the in- 
teraction of charged grains with the plasma becomes im- 
portant. The ion-grain cross section due to long-range 
Coulomb forces is larger than the atom-grain cross sec- 
tion. As a result, the rate of translational motion damp- 
ing gets modified. For subsonic motions the effective 
damping time due to gas drag is rcnormalizcd: 



a ^Td„ 



(5) 



with the following rcnormalizing factor (Draine 
Salpeter 1979) 



1 



On _ f ^ ^ \ TTi,^^ 

3 



Z 



\akBTt 



X In 



2\Z\^'Kxnii 



knT,, 



3/2 



(6) 



Here Xi is the abundance of ion i (relative to hydrogen) 
with mass and temperature Ti, x = '^^Xi, Ze is the 
grain charge, and fziZ) is the grain charge distribution 
function. When the grain velocity Vd relative to gas be- 
comes supersonic, the dust-plasma interaction is dimin- 
ished, and the damping rate in this case is renormalized 
due to the gas-dynamic correction (Purcell 1969), 



a = 1 



128 C2 



1/2 



(7) 



where Cs = \J k^Tnlmn is the sound speed. 

If tl is greater than r^ragi then the effect of magnetic 
field on dust dynamics is negligible. However, tl ^ Tdrag 
in most phases of the ISM. 

3. RESONANCE ACCELERATION 

3.1. Gyroresonance acceleration: nonlinear theory 

In this section, we revisit the treatment of resonance 
acceleration by fast modes in compressible MHD turbu- 
, (4)ence using nonlinear theory (NLT). 
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TABLE 1 

Idealized Environments and MHD turbulence parameters 



PdTCLTTlGtGTS 


CNM 


VV IN iVi 


WIM 


nil (cm^'^) 


30 


0.4 


0.1 


^gas (K) 


100 


6000 


8000 


XH 


0.0012 


0.1 


0.99 


B ii^G) 


6 


5.8 


3.35 


L (pc) 


0.64 


100 


100 


SV = Va (km s-i) 


2 


20 


20 


kc (cm~l) 


7 X 10^15 


4 X lO-i'^ 




Damping 


Neutral-ion 


Neutral-ion 


Ion \ 



and coUisionless 



"Here rtfj is the gas density, Tgas is the gas temperature, xjj is the 
ionization fraction of H, _B is the strength of magnetic field, L is the 
injection scale of turbulence, &V is the rms velocity of turbulence 
at the injection scale, and kc is the cutoff scale of turbulence due 
to coUisional and coUisionless damping. 



Consider a grain of mass m, charge Ze, moving with 
velocity v m. a, magnetized turbulent medium with a uni- 
form magnetic field B = Bq. The motion of such charged 
grain in B consists of the gyration of the grain about its 
guiding center and the translation of the guiding cen- 
ter along B. In the QLT limit, the guiding center is 
assumed to follow a regular trajectory along B with a 
constant cosine of pitch angle /i = cos /3 with /3 being 
the angle between v and B. Gyroresonant interactions 
between grain and wave occur when the wave frequency 
in a reference system fixed to the grain guiding center is 
a multiple of the Larmor frequency: 



fc||W/i 



ill 



(8) 



Tg or Z < Tg, where the 



with n = ±l,±2,.... This resonance condition is equiva- 
lently described by a Delta function 5„(a; — — n^l). 

Gyroresonance accelerates grains in the direction per- 
pendicular to the mean magnetic field Bq because electric 
field induced by plasma perturbations is perpendicular 
to Bq (see e.g. YLD04). This acceleration mechanism 
is dominant by eddies smaller than the grain gyroradius, 
i.e., I < Tg. Indeed, consider gyroresonance by fast MHD 
modes in low-/? plasma. From the resonance condition 
(8) for n = 1 with lu = fc||WA, we obtain turbulent scales 
for gyroresonance k > k^c-. 
fact that /i > — 1 has been used 

In a turbulent medium, the local magnetic field B = 
Bq + (5B where S'B is the turbulent component of mag- 
netic field, varies both in space and time, so fi changes, 
and U|| and v± change accordingly. The grain guiding 
center has fluctuations from its regular trajectory along 
Bq. The NLT takes into account such fluctuations of the 
guiding center. 

Assuming that the projection of the fluctuations of 
grain guiding center onto the mean field Bq can be de- 
scribed by a Gaussian distribution, the resonance condi- 
tion becomes 
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where n = and ±1, and Ad|| is the dispersion of velocity 
(see YL08 and YLP08; Appendix C). 

We arc interested in the grain acceleration, so the dif- 
fusion coefficient arising from gyro-phase averaging Dpp 
is used (sec Schlickeiser & Miller 1998). In compress- 
ible MHD turbulence, the fast modes are shown to be 
dominant in gyroresonance acceleration (YL03). Its cor- 
responding diffusion coefficient is given by (see Appendix 
C) 
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where n = ±1. In the above equation, L is the injection 
scale of turbulence, w = k±^v^/VL, x = fc/fcmin = kL, R = 
wfcniin/f2, Af^ = 6V^/Vj^, rj = cosO with 9 is the angle 
between the wave vector k and the mean magnetic field, 
kc is the cut-off of turbulence cascade due to damping, 
and J„ is second order Bessel function. The dispersion 
of the cosine of pitch angle A/i is given by Equation (C5) 
in Appendix C. 

3.2. Transit-Time Damping (TTD) 

Transit-time damping (TTD) or transit-time accelera- 
tion, arises from resonant interactions of particles with 
the compressive component of magnetic fiuctuations, i.e., 
the component parallel to the mean magnetic field Bq 
in magnetized turbulent environments. When the grain 
moves together with the wave along Bq, it is subject 
to magnetic mirror forces — (mi'^/2_B)V||B, where v±^ is 
the grain velocity component perpendicular to Bq. In 
the plasma reference, the back and forth collisions of the 
grain with the moving magnetic mirrors increase grain 
energy because the head-on collisions are more frequent 
than trailing collisions (see e.g., Fisk 1976). The reso- 
nance condition for a grain with velocity v reads 



to ■ 



0, 



(11) 



where v\\ = vfi is the grain velocity component paral- 
lel to B, and u! is the wave frequency (see Fisk 1976; 
Schlickeiser & Miller 1998). 

For fast MHD modes in low-/3 plasma, the dispersion 
relation is cj = kVA (see Cho et al. 2002), and the re- 
quired velocity for TTD corresponds to W|| = va/cosO. 
Thus, if W|| > Va, TTD can be efficient to accelerate 
grains to large velocities. 

In the NLT limit, the diffusion coefficient for TTD is 
given by (see Appendix C) 



TTD 
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drj 



(A/i)2 



(12) 



In Figure 3 we present Dp™ as a function of the cosine 
of the grain pitch angle fi obtained using the QLT and 
NLT. Two values of the grain velocity v = 1.5 Va and 
3.5 Va are considered. It can be seen that in the former 
(9)case, D'IJ^ increases with decreasing ^ until fi = v/Va, 



graphite, WIM 



n ' ' ' I ' ' ' r 
v=1.5V, 
v=3.5Va 




Fig. 3. — Diffusion coefficient Djp^ as a function of tlie cosine 
of tlie grain pitcli angle /i for the limit of QLT (green lines) and 
NLT (red lines). Two velocity values of super- Alfvcnic graphite 
grains of size a = 1Q~^ cm are considered. Djp^ drops sharply 
for fi <Va/v in the QLT, but Dj™ is finite as -i> in the NLT 
as a result of broadening of resonance conditions. 

and drops sharply to zero for fi < v/Va because the res- 
onance condition (11) is not satisfied. In contrast, the 
broadening of the resonance condition in the latter case 
allows grains with fj, < Va/v to have resonant interac- 
tions with waves, resulting in finite Dj™ even at ^ = 0. 
Therefore, TTD is usually neglected in the QLT because 
the gyroresonance tends to accelerate grains in the per- 
pendicular direction to Bq, resulting in /i = 0, for which 
D™ 0. However, TTD can play an important role 
in driving grain motion when the fluctuations of guiding 
center are taken into account. 



4. GRAIN VELOCITIES DUE TO RESONANCE 
ACCELERATION 

4.1. Grain dynamics 

Consider an ensemble of grains with the same mass m, 
moving in the uniform magnetic field with their different 
pitch angles ^. From the equation of motion, mdv/dt 



-mv/tdr 

obtain 



R, where R is the random force, we can 



dt 



^dra 



A{v), 



(13) 



where (w^) is the grain velocity dispersion averaged over 
the ensemble of grains, A{v) is the rate of energy gain 
(see YL03). 

When the scattering is less efficient than the accel- 
eration,^ the cosine of the grain pitch angle /i = cos (3 
changes slowly during acceleration, and A{v) is given by 



Aiv) 
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where Dpp(j), ^) is the diffusion coefficient. 



^ The efficiency of scattering relative to acceleration is de- 
scribed by the ratio p^-D^^/i)pp = cos {v cos 9/Va + Z^)^ where 
cos 6 = fcy/A; with 6 is the angle between k and B (see Appendix 
D). The scattering is negligible for grains with v < Va, but it be- 
comes important for v > Va or when the grain is moving along 
the magnetic field, i.e., tiii 2> v±- 



E(k) 




Fig. 4. — A sketch of spectrum of MHD turbulence and accel- 
eration mechanisms for charged grains are shown with its corre- 
sponding scale, fcinin ~ is the injection scale, l^^fp is the grain 
mean free path, fcros ^ ^'g^ is the gyroresonance scale, and is the 
damping cut-off of turbulence. Wc assume that grains of interest 
arc large enough so that fcres < fcc or Larmor period larger than 
eddy turnover time. 

When the scattering is more efficient than the accelera- 
tion, the pitch angle can be rapidly redistributed through 
pitch angle diffusion during the acceleration (i.e., diffu- 
sion approximation). The rate of energy gain (Eq. 14) is 
then determined by the averaged value Dp of Dpp{p,^) 
over the isotropic distribution of /.i (see Dung & Schlick- 
eiser 1990ab): 



-Dp(p) = ^ j ^ (^Dpp{p,fi) 



Dt,t,ip,tJ-) 



d^i. (15) 



We consider in the present paper zero helicity turbulence 
with Df^p = 0. When Dpp is known, we solve Equation 
(13) iteratively to get convergent velocities. 

Physical parameters for ISM conditions are shown in 
Tabic 1. MHD turbulence is injected at a large outer 
scale L with velocity dispersion 5V . The value 5V is cho- 
sen such that the turbulence is weak and sub-Alfvcnic. 
Major damping processes are also listed in Table 1. For 
the CNM and WNM, the dominant damping arises from 
the neutral-ion viscosity. The value of fee for these phases 
is adopted from YLD04. For the WIM, the turbulent 
damping arises mainly from ion viscosity and collision- 
less damping. We calculate by equating the damping 
rate to the rate of turbulence cascade (see Appendix C, 
also YL03; YLD04). 

Figure 4 sketches possible acceleration mechanisms 
present in MHD turbulence. Working scales in the in- 
ertial range, spanning from the injection scale fcmin to 
the cutoff scale kc of turbulence, are indicated. Gyrores- 
onance works only if fej-cs < kc, and is dominant at small 
scales k > fcj-cs ^ The critical size Ocri for gyroreso- 
nance is then obtained by solving the equation k^cs = kc 
for grain size. 

Figure 5 shows the rate of energy gain (Eq. 14) as a 
function of grain velocity v for gyroresonance accelera- 
tion (n = ±1) and TTD acceleration (n = 0) in the CNM 
{upper panel) and WIM {lower panel) arising from fast 
modes in MHD turbulence. Here we assumed that the 
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Fig. 5. — Rate of energy gain as a function of the grain velocity 
for gyroresonance acceleration (n = ±1) and transit time damping 
acceleration (n = 0, TTD) for a graphite grain of size a = 10~^ cm 
in the CNM {upper) and WIM (lower). Solid and dot lines denote 
results from NLT and QLT, respectively. The Alfvenic speed Va 
is indicated. The gyroresonance acceleration (n = 1) is dominant 
for 11 < Va, and the TTD acceleration (n = 0) takes over when 
V > Va- The case of efficient pitch angle scattering is considered. 

pitch angle scattering is efficient. Due to the broadening 
of resonance condition, gyroresonant interactions occur 
at lower velocities in the NLT than QLT (solid lines). 
It can also be seen that the gyroresonance acceleration 
(n = 1) is dominant for v < Va, while TTD acceleration 
(n = 0) becomes dominant for v > Va- 

In both the CNM and WIM, as v ^ Va, = 1) 

becomes smaller in the NLT than the QLT. As a result, 
we expect that gyroresonance acceleration is less efficient 
in the former case. We are going to quantify such a 
difference in the following. 

4.2. Gyroresonance in quasi-linear theory and nonlinear 

theory 

Grain velocities due to gyroresonance acceleration are 
obtained by solving Equation (13) using the diffusion co- 
efficients from Equations (10) and (15). We assume that 
at the beginning the grain has low velocity, so that the 
pitch angle scattering by TTD is negligible. The gyrores- 
onance increases rapidly v±, and fi decreases to /i = 0. 
So, we can assume /i = for the gyroresonance acceler- 
ation. 

Figure 6 shows grain velocities obtained using the NLT 
and QLT for the CNM, WNM and WIM. Both silicate 
and graphite grains are considered. Grain velocities ob- 
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Fig. 6. — Grain velocities relative to gas as a function of grain 
size for graphite and silicate grains in various phases of the ISM. 
Acceleration arising from gyroresonant interactions of fast modes 
with grains are obtained using the QLT (dotted line) and NLT 
(solid line). The difference in velocities from the NLT and QLT is 
within 15%. The shaded area indicates the range of size in which 
the assumption of constant charge is invalid due to strong charge 
fluctuations. Dotted vertical lines present in the WNM and CNM 
denote the critical size Ocri corresponding to the cutofi scale fcc of 
turbulence. 

taincd using the NLT are generically smaller than those 
from the QLT. But the difference is within 15%. The 
smaller results in the NLT arise from the fact that, when 
grain velocities approach V^a, a fraction of turbulence 
energy is spent to induce transit time acceleration (see 
also Fig. 5). The sudden cutoffs (dotted lines) present 
in the CNM and WNM correspond to the cutoff scales 
of turbulence due to coUisional and collisionless damping 
occurring at the critical size (see Sec. 4.1). 

Figures 7 and 8 compare grain velocities arising from 
fast MHD modes (similar data as Fig. 6) with those 
induced by Alfvcn hydro-drag modes (LY02) and fast 
hydro-drag modes (YLD04) in the CNM and WIM, re- 
spectively. ^As expected from earlier studies, gyrores- 
onance acceleration is dominant for the entire range of 
grain size in the WIM (Fig. VNQwim). The rapid vari- 
ation of grain mean charge (Z) present in the WIM (see 
Fig. 1) results in the non monotonic increase of grain ve- 
locities from gyroresonance (solid and dotted lines). Two 
local maxima in the hydro drag cases correspond to the 
change in sign of (Z) (see Fig. 1). In the CNM, gyroreso- 

^ Here only the velocity component perpendicular to magnetic 
held is shown. 
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Fig. 7. — Comparison of grain velocities arising from gyroreso- 
nant interactions of fast MHD modes {solid and dotted lines) with 
the results arising from hydrodrag by fast and Alfven modes {dot- 
dashed and dashed lines). 

nance is dominant for grain size from Ocri = 5 x 10^® cm 
to a ~ 6 X 10~^ cm, while acceleration by hydro-drag 
takes over for grains smaller than Ccri and larger than 
~ 6 X 10-5 cm (sec Fig. 7). 

4.3. Acceleration by TTD 

Gyroresonant acceleration tends to drive grain motion 
in perpendicular direction to the mean magnetic field 
(i.e., /i = 0). As discussed in Section 3.2, the broaden- 
ing of resonance condition in the NLT allows TTD to 
operate even at /i = 0. In addition to acceleration, the 
scattering by TTD can be important, which results in the 
deviation of ^ from fi = (see Yan & Lazarian 2008). 
However, it is still uncertain how fast the pitch angle 
scattering (both gyroresonance and TTD) by fast MHD 
modes is compared to the acceleration. For simplicity, 
we consider the TTD acceleration for two limiting cases 
of efficient scattering and inefficient scattering in which 
the scattering is more and less efficient than the acceler- 
ation, respectively. In the latter case, the scattering is 
assumed to be sufficient to alter the adiabatic invariant 
of gyromotion. 

In the presence of TTD, we take into account the diffu- 
sion coefficients D^J^ from Equation (12) for Equation 

(13) in addition to D^^. 

When the pitch angle scattering is efficient, ji is de- 
scribed by an isotropic distribution f{ii)diJL — l/2d^. 
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Fig. 8. — Similar to Fig. 7 but for the WIM. The shaded area in- 
dicates the range of grain size in which the assumption of constant 
charge is invalid. Local maximum velocity for hydro acceleration 
present in the range a < 10^^ and a > 10^^ cm arise from the 
change in sign of {Z). 



The diffusion coefficient Dp(j),^) is replaced by its av- 
erage value over the isotropic distribution. f{^). 

Figure 9 compares grain velocities relative to gas ob- 
tained using the NLT from gyroresonance {dot line), and 
gyroresonance plus TTD by fast MHD modes for silicate 
grains (upper) and graphite grains (lower) assuming the 
efficient pitch angle scattering. In the WIM conditions, 
the TTD acceleration is negligible for grains smaller than 
- 4 X lO-^cm for which v <Va- The efficiency of TTD 
begins to increase rapidly with a when v ^ Va- TTD can 
increase grain velocities to an order of magnitude higher 
than gyroresonance. The effect of TTD is less important 
in the WNM than the WIM, but still considerable. For 
the CNM, TTD acceleration is rather marginal because 
grains are moving with sub-Alfvenic velocities, less than 
the threshold for TTD. 

In the case of inefficient pitch angle scattering, the co- 
sine of the grain pitch angle /i = is assumed as a result 
of gyroresonance acceleration. We found that grain ve- 
locities are ~ 10% larger than the results for the case of 
efficient scattering. This seems counterintuitive because 
the scattering is required to alter the adiabatic invariant 
in gyromotion. However, here we neglected that effect of 
pitch angle scattering, and the situation is merely related 
to the fact that the diffusion coefficient Dpp averaged over 
fi is slightly lower than Dpp at /x = (see Fig. 3). 
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Fig. 9. — Grain velocities relative to gas arising from gyroreso- 
nant acceleration {dotted line) and gyroresonant acceleration plus 
transit time acceleration (solid line), as a function of grain size for 
different environments. The case of efficient scattering is consid- 
ered. 

5. STOCHASTIC ACCELERATION BY LOW FREQUENCY 
ALFVEN WAVES 

5.1. General consideration 

The stochastic acceleration by low frequency Alfven 
waves with uj < at the gyro-scale k± ^ p^^ was stud- 
ied in Chandran et al. (2010) for ion heating in the solar 
wind. Here we consider the effect of low frequency Alfven 
waves on dust acceleration in the ISM. For the low fre- 
quency Alfven waves, resonance acceleration is inefficient 
because the resonance condition w — k\\v^ = nfi is not 
satisfied. Indeed, in low-/3 plasma, the ion thermal veloc- 
ity vt ^ V^A, and lj ~ fc||V^A for Alfven waves, we have 
to — k\\v^ <C 0. 

Let 5v and 5B be the rms amplitudes of velocity and 
magnetic field at the gyro-scale k±p ~ 1. The electric 
field induced by plasma perturbations with velocity 5v in 
the direction perpendicular to the mean magnetic field 
has the rms amplitude 



5E'- 



SvBo 



(16) 



and the potential is written as 

6^ ~ pSE. 



(17) 

The electric field in the plane perpendicular to Bq re- 
sults in the acceleration of grains in the perpendicular 
direction. 



Combining these above equations, we obtain 
p6vBo 

qd<P ~ q = mv±ov, 



(18) 



where m is mass of the charged particle and p = 
mcv±/qBQ is the gyro radius. 

The fractional increase of energy after a single gyro- 
period is calculated by 



qS^ 
mv\l2 



2e, 



where 



(19) 



(20) 



When e ^ 1, also corresponding to 8B <^ Bq due to 
the assumption of low /3 plasma, the increase of energy 
per gyroperiod is negligible because of adiabatic invariant 
for the magnetic moment mv^/2Bo. As e increases to 
unity, the particle energy changes substantially after one 
gyroperiod. The guiding center becomes chaotic when e 
exceeds some threshold value, so that the perpendicular 
acceleration becomes important. 

The increase of energy depends only on the amplitude 
of perturbation at the scale of gyroradius. The increase 
of energy per unit of time per unit of mass is defined as 



Qa 



'Tacc 



(21) 



where Tacc is the time for the particle kinetic energy K± 
to increases by a factor of 2. We can estimate Tacc using 
the diffusion coefficient as follows: 
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Dk 



2„4 



m V 



D 



K 



where 



At 



2 2 2 r 



(22) 



(23) 



and Weft = Sv/p is the effective frequency of gyroscale 
fluctuations. Here we have taken AK± ~ mv±6v = 
mv^ujcsp from Equation (18), and At = p/Sv is the 
time required for the guiding center to move by a dis- 
tance equal to the gyro radius p. 

Taking use of Equations (22) and (23) for (21), we 
obtain 



Q± - ^cS^VP - 



(24) 



When e is sufficiently small, the variation in the parti- 
cle energy is correlated over long time, so the assumption 
for diffusive approximation may not be adequate. Thus, 
the actual energy gain is substantially smaller than the 
value obtained by Equation (24). To account for this re- 
duction, a damping function cxp(— C2/e) is introduced, 
and Equation (24) can be rewritten as 



Qj 



p 



-exp 



(25) 



where ci and C2 arc dimcnsionless constants that depend 
on the nature of fluctuations. Below we assume ci = 
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0.75 and C2 = 0.34 for the fluctuations by low frequency 
Alfven waves as in Chandran et al. (2010). 

5.2. Grain velocities for the ISM 

Let assume that the Alfvenic turbulence in the ISM 
follows the scaling 



5v 



(26) 



where Z_l ~ pd and L are gyro scale and the injection 
scale, a and a = (cs — l)/2 with C3 being the slope 
of turbulence power spectrum are dimensionless, derived 
from the properties of turbulence. For sub- Alfvenic tur- 
bulence, a < 1. Using Equation (20) for dust grains, i.e., 
ea = 5v_i_/v^_cU we obtain 



(l-a)/2 ^(l+a)/2 



(27) 



where p denotes proton, and d denotes dust, A = m/nip, 
dp = Va/^^p, and the perpendicular temperature is de- 



fined as kT±_ 



Using Equations (26) and (27) for Equation (25), we 
obtain the rate of energy gain per a grain of mass m 



T^exp 



(28) 



For calculations, we assume C3 = 5/3 and a = 1/3 for 
Alfvenic turbulence above the gyro scale. 

Using the parameters for the ISM in Table 1, we cal- 
culate the grain velocity arising from the chaotic accel- 
eration by low frequency Alfven waves in Figure 10 for 
the CNM and WIM. We show that the chaotic acceler- 
ation by low frequency Alfven waves is subdominant to 
the fast and Alfvenic hydrodynamic drag. Obviously, it 
is much less important than gyroresonance and TTD by 
fast modes. The possible reason is that the low frequency 
Alfven waves cascade faster to small scale than the fast 
modes. 

6. DISCUSSION 
6.1. Related works on dust grain acceleration 

The acceleration of dust grains by incompressible MHD 
turbulence was first studied by Lazarian & Yan (2002). 
Yan & Lazarian (2003) studied grain acceleration in com- 
pressible MHD turbulence, and discovered a new accel- 
eration mechanism based on gyroresonant interactions 
of grains with waves. This acceleration mechanism in- 
creases grain velocities in perpendicular direction to the 
mean magnetic field. YLD04 computed grain veloci- 
ties arising from gyroresonance by fast MHD modes us- 
ing quasi-linear theory (QLT), and compared the ob- 
tained results with different mechanisms, for various ISM 
phases. They found that the gyroresonance is the most 
efficient mechanism for grain acceleration in the ISM. 

The effect of large scale compression on grain accelera- 
tion is shown by Yan (2009) to be less important than the 
gyroresonance in the ISM conditions, unless the grains 
move with super-Alfvenic velocities. 

For very small grains (e.g., polycyclic aromatic hydro- 
carbons and nanoparticles), Ivlev et al. (2010) sketched 
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Fig. 10. — Grain velocity due to stochastic acceleration by low 
frequency Alfven waves (solid lines) compared to the acceleration 
by hydrodynamic drag from Alfven modes (dotted lines) and fast 
modes (dashed lines), for the CNM {upper) and WIM {lower). The 
stochastic acceleration is much less efficient than the latter. 

a new mechanism of grain acceleration due to electro- 
static interactions of grains with fiuctuating charge and 
provided rough estimates of grain velocities in the ISM. 
Hoang & Lazarian (2011) quantified this mechanism us- 
ing Monte Carlo simulations of charge fluctuations. They 
found that charge fluctuations can accelerate grains to 
several times their thermal velocities. 

6.2. NLT for gyroresonance acceleration 

We have revisited the treatment of gyroresonance ac- 
celeration for charged grains due to MHD turbulence by 
accounting for the fluctuations of grain guiding center 
from a regular trajectory along the mean magnetic field 
(i.e. NLT limit). The fluctuations of the guiding center 
result in the broadening of resonance conditions- a Delta 
function is replaced by a Gaussian function. Such broad- 
ening of resonance condition allows some fraction of wave 
energy spent through the TTD acceleration. As a result, 
grain velocities due to gyroresonance acceleration are in 
general decreased by ~ 15% in the NLT limit. 

6.3. Transit time damping acceleration 

TTD acceleration is believed to be important when the 
parallel component of grain velocity along the magnetic 
field exceeds the Alfven speed Va- Although gyroreso- 
nance acceleration by fast modes can accelerate grains to 
V > Va, their resulting velocity mostly perpendicular to 
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the magnetic field, i.e. /i = 0, makes TTD mifavored be- 
cause the resonance condition (5(aj — fc||W||) is not satisfied. 
Indeed, we found that TTD is efficient for fi > Va/v and 
negligible for /i < Va/v in the QLT limit. This feature is 
consistent with the result for acceleration of cosmic rays 
in Schhckciser & Miller (1998). 

The situation changes when the fluctuations of the 
guiding center are taken into account in the NLT. For 
this case, the resonance condition is broadened beyond 
the S function, and can be described by a Gaussian func- 
tion. As a result, TTD acceleration becomes important 
for /i < Va/v, including 90° pitch angle. 

In addition to acceleration, TTD also induces the grain 
pitch angle scattering, which is dominant over the scat- 
tering by gyroresonance. Since the efficiency of the TTD 
scattering is uncertain, we considered in the paper two 
limiting cases of inefficient and efficient scattering in 
which the scattering is less and more efficient than the 
acceleration. The pitch angle is equal to 90° in the for- 
mer, and isotropic in the latter. 

When the scattering is more efficient than the acceler- 
ation, we showed that for the WNM and WIM, the TTD 
acceleration can increase substantially the grain velocity 
compared to results arising from gyroresonance. Particu- 
larly, for grains larger than 5 x 10~^ cm in the WIM, TTD 
acceleration is an order of magnitude greater than the gy- 
roresonance acceleration. TTD is clearly more efficient 
than the betatron acceleration studied in Yan (2009). In 
the CNM, TTD acceleration is limited because the gy- 
roresonance acceleration is not able to speed up grains 
to super Alfvenic stage. When the scattering is less ef- 
ficient than the acceleration, grain velocities are within 
10% lower than the results for the efficient scattering 
case. 

6.4. Stochastic acceleration 

We study also the effect of low frequency {uj < ^l) 
Alfvcn waves on dust grains in the ISM. We show that 
the stochastic acceleration by low frequency Alfven waves 
is subdominant to the gyro-resonance acceleration and 
TTD acceleration by fast modes. This may arise from 
the fact that low frequency Alfvcn waves cascade faster 
to small scale than the fast modes in MHD turbulence. 

6.5. Implication to dust coagulation and shattering and 

alignment 

Hirashita & Yan (2009) adopted grain velocity due to 
gyroresonance acceleration from YLD04 to model grain 
size distribution in the different ISM conditions. Hi- 
rashita et al. (2010) studied grain coagulation and shat- 
tering in the WIM for large dust grains ejected from Type 
II supernova. They showed that the shattering of large 
dust grains due to turbulence plays an important role in 
producing small size population that modifies extinction 
curves in starburst galaxies. 

The threshold velocity for the grain shattering is a 
function of the grain size: 

Wshat = 2.7 (— -^^ ) "^\ms-\ (29) 

\1G~' cm/ 

where ?;dd is the relative velocity of dust grains (Chok- 
shi et al. 1993). If Udd < Vshat, the grains collide and 
stick together. When Vdd > ^shat; the collisions with 



high velocity produce shock waves inside the grains, and 
shatter them in smaller fragments. For Vdd > 20 km/s, 
the evaporation of the dust grain occurs and the grains 
are destroyed. 

When TTD is accounted for, grains larger than 5x10^^ 
in the WIM and WNM may undergo efficient shattering 
because vm > V^a = 20 km s~^. 

The effects of high velocities we obtain on dust grain 
alignment require further studies. Recent research has 
shown that the classical grain alignment theory of Davis 
& Greenstein (1951) (see also Lazarian 1995; Roberge 
& Lazarian 1999 for more recent quantitative studies 
of the process) is subdominant to the radiative torque 
(RAT) model (see Dolginov & Mitrofanov 1976; Draine 
& Weingartner 1996; Lazarian & Hoang 2007a). This 
model, however, was criticized in Jordan & Weingartner 
(2009) who appealed to the results of gyroresonance ac- 
celeration of grains in Yan & Lazarian (2003) and claimed 
this means that fast moving grains will be randomized as 
their charge fluctuates. Our results indicate that grains 
can be accelerated to even faster velocities, which could 
make the problem for the alignment more severe. How- 
ever, we believe that the claim about suppression of the 
RAT alignment for fast moving grains is a result of the 
confusion on the nature of the RAT alignment. We 
plan to address this issue elsewhere. At the same time 
high velocities of grains may induce another mechanical 
alignment of irregular grains as it described in Lazarian 
& Hoang (2007b). This alignment does require further 
studies. 

7. SUMMARY 

In the present paper, we study the resonance accel- 
eration of charged grains by fast modes and stochastic 
acceleration by Alfven waves in MHD turbulence. Our 
main results are summarized as follows. 

1 . We revisit the treatment of gyroresonance accelera- 
tion of charged grains in compressible MHD turbulence 
by taking into account the fluctuations of grain guiding 
center from the regular trajectory along the mean fleld. 
Gyroresonance interactions by fast modes can acceler- 
ate large grains to super- Alfvenic speed. We found that 
grain velocities are lower by 15% in the NLT than the 
QLT. 

2. We investigate the effect of transit time damp- 
ing (TTD) by fast modes for super- Alfvenic grains. We 
found that the fluctuations of grain guiding center allow 
TTD to occur not only within the range of the cosine of 
the grain pitch angle ^ > Va/v as expected by the QLT, 
but also for ^ < Va/v. We show that the TTD acceler- 
ation can increase grain velocities by an order of magni- 
tude compared to the results arising from gyroresonance 
mechanism Thus, TTD is the most efficient acceleration 
mechanism for super- Alfvenic grains. 

3. The stochastic acceleration due to low frequency 
Alfvcn waves is inefficient for dust grains in the ISM con- 
ditions. 
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APPENDIX 

A. TURBULENCE CASCADE AND DAMPING 

We summarize here major damping processes for the MHD turbulence. First, we begin with the cascade of different 
MHD turbulence modes. 

A. l. Turbulence cascade 

Turbulence cascades from large scale to small scale. The cascade rate depends on the scale of the eddy and on the 
type of turbulence mode. For Alfven and pseudo-Alfven (slow) modes, the cascade time at the scale k is equal to the 
eddy turn-over time: 

^^"^^ = ^kl^- ^^^^ 

where the critical balance condition (GS95) for the cascade parallel and perpendicular to the mean magnetic field 

A:||Wa = k±v± (A2) 

has been used. Here v±_ ~ Vk because the Alfven mode has perturbation velocity Vk perpendicular to the mean 
magnetic field. 

Fast modes cascade a bit slower than the eddy turn over, and the cascade time is given by 

where Vph is the phase speed. In low-/3 plasma, Vph = v a and Vph = Cg in high-/3 plasma. 

B. TURBULENCE DAMPING 

B. l. Collisional damping 

B.1.1. Neutral Ion collision damping 

In partially ionized gas, the MHD turbulence energy is dissipated through ion-neutral collisions. Let A„ be the mean 
free path of a neutral. For scales / > A„, the ion- neutral collisions can damp the turbulence at a damping rate 

r,„; Vnl^"^ ~ ^^(A„w„)fc^, (Bl) 

^gas 

where / ^ fc^^ is the scale of interest, Vn is the effective viscosity due to neutrals, and rigas — Ui + Un is the gas number 
density. 

B.l. 2. Viscous damping 

In fully ionized gas, the turbulence can get damped due to viscous damping or coUisionless damping depending on 
whether regimes are collisional or coUisionless. The Coulomb mean free path of the thermal plasma is given by 

V m^v-^ 7 / Tgas \ VlO^°cm-n 
Acoui = ^^Tcoui = = T = 9 X lOM cm, B2 

where cr = Trrf^^^^ = tt (e^/mw^) has been used. 

The viscous damping is important for scales between the injection scale L and the mean free path Acoui- For scales 
I < Acoui, the coUisionless damping is dominant. 

B.2. CoUisionless regime 
B.2.1. ion viscosity damping 

The motion in the perpendicular direction to the magnetic field is suppressed, so the viscosity is much smaller 
in the perpendicular direction compared to the viscosity in the parallel direction, i.e. ryj_ ~ ?7o/(^^iTCoui)^ where 
T'Coui ^ ^Coui/vth is the Coulomb collision time for ion and Vth is the ion thermal velocity, and r]o = O.QGnfcsTrcoui is 
the longitudinal viscosity for motion along the magnetic field (see Braginskii 1965). Petrosian et al. (2006) derived 

brumi 

for /3 < 1 and 
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for /3 ^ 1. For fast modes, the cut-off due to viscous damping is obtained in YL08: 
for /3 ^ 1 and 

kcL = Xcil - 3 cos ^^)"'^/^ 

for the high /3 plasma, where Xc = (6pSV^L/{r]QVA))^^'^ ■ 
For fast modes, the cut-off scale of turbulence is given by 



Lf3 



-1/3 / \ ■ n2 \ -2/3 

Ml 



4 /A/aV^V L 
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(B5) 
(B6) 

(B7) 



B.2.2. Landau damping 

In collisionless plasma, the turbulence is damped due to the Landau damping. The damping rate due to the Landau 
damping is given by (see YL08) 

1/2 
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(B8) 



The cut-off scale of the turbulence is obtained by equating the damping rate to the cascading rate. From Equations 
(A3) and (B8) we obtain 
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(B9) 



7rTOe/3sin( 

C. NLT OF GYRORESONANCE ACCELERATION OF DUST GRAINS 
Below we describe the NLT for resonance acceleration of dust grains in MHD turbulence. 

C.l. Fokker-Plank coefficients 

The QLT assumes that the guiding center of the charged particles is regular motion along the uniform magnetic field 
and that the gyro-orbit is not perturbed. 
The Fokkcr-Planck diffusion coefficients (Jokipii 1966; Schlickeiser & Miller 1998) are given as 



pp 
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where ub = BI/%tt, \kmin\ = K 
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(CI) 



corresponds to the dissipation scale, TZ, C refer to the left- 



and right-circularly polarized modes, and (j) ~ tan^^ k^/ky. Above i?„ is the function for resonance condition, Vph is 
the phase speed and D, is Larmor frequency. 
The correlation tensors are defined as 



(B„(k,t)m(k',t + T))/B2 =5(k-k')A/«/3(k). 



t)v^k\ t + t))IVI = S{k - k')A-,Mk)e 



j-r/Tfc 
--r/Tfc 



(C2) 



where Ba,p, Va,p are respectively the magnetic and velocity perturbation associated with the turbulence, is the 
nonlinear decorrelation time and essentially the cascading time of the turbulence. For the balanced cascade we consider 
(see discussion of our imbalanced cascade in CLV02), i.e., equal intensity of forward and backward waves, Cij(k) = 0. 
The magnetic correlation tensor for Alfvenic turbulence is (CLV02), 



M,,(k) 
i^„(k) 



127r 



/,fc;^"/^xp(-LV3|^.|il/fcf), 



(C3) 



where lij = {Sij — kikj/k'^} is a 2D tensor in x ~ y plane which is perpendicular to the magnetic field, L is the 
injection scale, V is the velocity at the injection scale. Slow modes are passive and similar to Alfvcn modes. The 
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normalization constant is obtained by assuming equipartition £k — J dk^ ^'luBQ/Sn ^ Bq/Stt. The normalization 
for the following tensors below are obtained in the same way. 

According to CL02, fast modes are isotropic and have one dimensional energy spectrum E{k) oc fc^^/^. In low /? 
medium, the corresponding correlation is (YL03) 



My(k) 



Stt 



cos" 

1 



n = (k/L)-'/^ X Va/V, 



(C4) 



where 9 is the angle between k and B, iJ.y = kikj/k'j_ is also a 2D tensor in a; — y plane. The factor cos^ 9 represents the 
projection as magnetic perturbation is perpendicular to k. This tensor is different from that in Schlickeiser & Miller 
(1998). For isotropic turbulence, the tensor of the form cx; Ek{Sij — kikj/k'^) was obtained to satisfy the divergence free 
condition k • S'B = (see Schlickeiser 2002). Nevertheless, the fact that S'B in fast modes is in the k-B plane places 
another constraint on the tensor so that the term Sij doesn't exist. 

C.2. Diffusion coefficients in NLT 

The motion of a charged particle in a magnetic field B consists of the motion of the guiding center with respect to 
the magnetic field B and the motion of the particle about the guiding center. In the QLT limit, the guiding center is 
assumed to follow regular trajectory with constant pitch angle fi. In MHD turbulence, B varies with respect to space 
and time, so fi changes, and and v±^ change accordingly. 

In the NLT limit, the dispersion of the pitch angle due to magnetic field fiuctuations reads 

A^i = ^, (C5) 



and 



Av\\ _ {{B - Bo 



,2x1/4 



B, 
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^0 
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B2 
^0 



(C6) 



where Bq is the mean magnetic field (Volk 1975). The dispersion of parallel velocity Sv^^ is mainly induced by the 
fluctuations of the parallel magnetic field while 6B± is only second order effect. 

Since fi is constant in the QLT, the resonance condition is given by 5{k\\fiv — w + nQ). In the NLT, due to the 
fluctuations of fi (Eq. C5), the resonance condition is broadened and described by 



Rn (knvfi 



in) ^ 



kuAvi 



-exp 



(fejlij/i — u! + nfiy 



where n = and ±1 (see YL08; YLP08). 
Using Equations (C4) for (CI) combined with (C7), we obtain 



drj 



^l(At'it)^ 



[Joiw) + Jl{w)]cxp 



(C7) 



(^_^± 1 )^ 

^' rjv rjxK ' 



(C8) 



for the gyro-resonant acceleration n = ±1, and d?'k = l-nk^dkck] has been used for fast modes. In the above equation, 
L is the injection scale of turbulence, w = fc^w^/f2,x = k/k^^ii^ — kL,R ~ vkmin/^, M\ = dV'^/Vj^, rj = cos9, kc is 
the cut-off of turbulence cascade due to damping, and J„ is second order Bessel function. 
For the transit time acceleration (TTD), n = 0, we obtain 



TTD 



2Li?2 



-'/^dx 



drj 



(u;)exp 



Va\2 

TjV ' 



where L is the injection scale of turbulence, w = /c^u^/il, x = k/kmin = kL, R 
In the QLT, i?„ is replaced by S function. Equations (C8) and (C9) become 



vk„ 



A^2 
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(C9) 
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